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IN VISCOUS FLOW OVER AIRFOILS* 

R. E. Melnik 
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SUMMARY 

This paper describes a new theory developed by the author and his 
colleagues for analyzing viscous flows over airfoils at high Reynolds 
numbers* The theory Includes a complete treatment of viscous Interaction 
effects Induced by the curved wake behind the airfoil and accounts for 
normal pressure gradients across the boundary layer in the trailing edge 
region. A brief description of a computer code that was developed to 
solve the extended viscous interaction equations is given. Comparisons of 
the theoretical results with wind tunnel data for two rear loaded airfoils 
at supercritical conditions are presented. 

INTRODUCTION 

Two classes of theoretical methods are currently employed in the analysis 
of viscous flows over airfoils; one based on a direct numerical solution 
of the Reynolds-averaged equations of turbulent flow (e.g., refs. 1 and 2) 
and the other on an approximate boundary layer type formulation. Because 
of the magnitude of the couq>utational task direct solutions of the Reynolds 
equations have been restricted to relatively crude grids. Because of the 
long computing times and poor resolution, these methods are not yet suitable 
for practical applications. Most existing methods of analyzing the problem 
are based on the displacement surface concepts of conventional boundary layer 
theory (e.g., see ref. 3). Unfortunately these approaches themselves are 
not entirely satisfactory because they neglect certain strong interaction 
effects that occur near trailing edges and shock Impingement zones. In 
addition, most of the boundary layer formulations employed in the past were 
incomplete in that they ignored certain effects due to the wake. These 
effects are formally the same order of magnitude (as a function of the 
Reynolds number, R for R «) as the usual displacement effects on the 
airfoil surface ana should be Included Ip a completely consistent theory. 

There are two wake effects, one due to its thickness and the other to its 
curvature. The wake thickness effect leads to a semi- infinite equivalent 
displacement body that enters Into the determlnaton of the outer Invlscid 
solution. Although there are no conceptual difficulties in doing this, the 
wake thickness terms have usually been dropped or drastically approximated 
(e.g., ref. 3) in most previous studies in the Interest of simplifying the 
computations. The wake curvature effect arises from the turning of the low 
momentum flow in the wake along the curved streamlines behind the airfoil . 

*Thls work was partially supported by NASA under Contract NAS 1-12426. 
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thin shear layers on the airfoil, and in the wake and Strong interaction 
regions near the trailing edge and shock inplngenent points. The conventional 
boundary layer forsmlation leads to matching conditions coupling the invlscld 
flow to the thin shear layers outside the strong interaction regions. There 
are several equivalent statements of the matching conditions (see ref. 10) 
that follow directly from the fact that the flow in the shear layers is 
governed by the boundary layer equations. The matching conditions are 
Independent of the particular closure assumptions used to model the Reynolds 
stresses, and the various forms are all equivalent In the limit R -*■ ». In 
the present study we have used the surface source formulation because It was 
best suited for matching to the strong Interaction solutions aud because 
It eliminates the need to modify the airfoil geometry In the solution 
procedure. 


In 

viscous 

V 

AV 


this formulation the Invlscld flow equations are solved subject to 
boundary conditions on the airfoil and in wake given by: 



dp U 6 * 
e. e 


d S 

dp V 6 * 
e e 

d S 


on the airfoil surface 
(N - 0) 

on the wake streamline 

(N - 0) 


( 1 ) 

( 2 ) 


AU 


U (6*+e)ic on the wake streamline 

® (N - 0) 


(3) 


where V is velocity component normal to the airfoil surface, AV and AU are the 
jumps in the normal and tangential components of velocity across the wake 
streamline, Pg and Ug are the density and tangential velocity along the airfoil 
and wake, &* and 6 are the displacement and momentum thicknesses of the shear 
layers, S and N are curvilinear coordinates along the normal to the streamline 
defining the airfoil and wake centerline, and k is the curvature of the wake 
streamline. The three matching conditions can be formally derived from a 
standard asymptotic analysis in the limit R^ 


Since the wake curvature condition is not as familiar as the other 
conditions we siimmarize the steps in its derivation. The condition arises 
from the pressure variation induced across the wake by the nonuniform velocity 
profile in the wake. A typical velocity profile and pressure variation in the 
wake are sketched in figure 2. In an invlscld flow the pressure gradient 
3p/3N \jould be roughly constant across the wake and equal to 


(t^) “ - P 

InvlsciS 


U 


(4) 


In the viscous flow the pressure gradient is related to the actual velocity 
and density profiles in the wake U(N,S), p(N,S) by the normal momentum 
equation, as follows 


(5) 


ip. 

3N 


p(N,S) U^(N,S)k(S) 
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assumed that the streamlines In the wake are parallel, a 
justifiable assumption only in the part of the wake outside the trailing edge 
region. Subtraction of equations (4) and (5) followed by Ir.tegratlon wlth°* 
respect to N and the imposition of the condition that the pressure in the 
viscous flov approach the linear variation implied by equation (4) as 
N ± 00 leads to the relation ^ ^ 


AP = p" 


- Cj(S)k(S) 


( 6 ) 


where 


- P U ' 
e e 


(6*+6) 


<0 


(7) 


* ^ limiting values of the pressure in the outer inviscid flow 

on the upper and lower sides of the wake streamline, as indicated in figure 
. Equation (6) is similar to the boundary condition that arises in let 
ap aM ys s, except that there Cj is the jet momentum coefficient which is 
a positive constant that is given as part of the data. 


relaterbrAP^-‘*^®T^^'Jn^‘'^r’vJ^® pressure and velocity discontinuities are 
^ Combining this expression with Equation (6) leads 

an irrotational approximation to the outer inviscid flow and the wake 

curvature condition is Implemented by prescribing a Jump S velocUy potential 
r, across the wake streamline, where velocity potential. 


r = 


Ilmly^^ ■ condition used in the computation. 


( 8 ) 


~ = - u^(6*+e) 


18 

dS 
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^ere e is the angle the wake streamline makes with the chord of the airfoil. 

il of the n^erical solution, the right side of equation (9) 

is known and F can be evaluated by simple quadrature. 

TRAILING EDGE CORRECTIONS 

AIT source velocity, V, streamline curvature, ic, and the velocity iumo 

AU, appearing in the matching conditions all develop square root singlllrities 
and become unbounded at the traUing edges of lifting airfoils (see ref. 7). 

i ’ I condition cannot be satisfied, the lift cannot be 

determined and the conventional boundary layer theory cannot strictly be used 
to determine the viscous flow over lifting airfoils. Existlnrlelhods 
employing boundary layer theory apparently circumvent these shortlomings by 
numerically smoothing the inviscid and boundary layer solutions near traillne 

'l-e coupled laulecld end bolder, UyreJ^tiJJe lo * 
Obtain a self-consistent solution. The iteration seams to eliminate the 

consistent with an order of magnitude analysis 
of the normal momentum equation. ^ 
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A completely consistent treatment of the flow near trailing edges was 

associates for airfoils with cuspfd trailing 
g 8 (see ref. 7). In that analysis we showed that normal pressure 

fludirr boundary layer and wake are important and must be in- 

ff description of the local solution. We also showed 

tational flow, and we obtained complete analytic solutions of the relevant 
equations governing the leading terms of the local solution. Readers 
detalH^^*^ aspect of the problem should consult ref. 7 for further 

In this section we indicate how the local trailing edge solution can 
be used to correct conventional boundary layer theory to obtain composite 
solutions that are uniformly valid near trailing edges. In our apwoach we 
represent the solution by standard composite expressions that are valid 
in the invlscid, boundary layer, wake and trailing edge regions. For 

-efficient Cp and normal velocity 


from the invlscid solution 


ACv 


S = Vouter^N;R^)V ICp^^^^^^(S,N;R^) - Cp^^^^(S,N;R^)] 
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trailing edge. The last term with subscript "com" is the common part of 
he outer and Inner solution and is determined as part of the inner solution 
We should stress that the individual terms in the representation are not 
expanded in formal sum type asymptotic series in terms of the Reynolds number 
8 in the usual method of matched asymptotic expansions. Instead, the inner 

Z sLaJ U^r eq'uat Iterative solution of the coupled invlscid 
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layer and invlscid solutions outside the trailing edge region and to the 

strong interaction solution near the trailing edge. The fact 
that this representation can be made to work is due to the close relation 
between the downwash generated by the vortlcity in the shear layers and the ^ 
surface source velocities generated by the displacement thickness. 
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The appropriate correction to the matching conditions are generated as 
follows: First, the source velocities, V, on the upper and lower sides of 

the airfoil are linearly combined Into symmetric and antisymmetric components 
defined by 

V 2 = 2 <12) 

where V and V are computed from equation (1) evaluated on the upper and 
lower sides of the airfoil, respectively. Then V and V_ are multiplied 

by corrections G. and AS 

A S 

corrected ^A ^A ’ corrected ~ 

where G and Gg are the ratios of the Inner solution to Its common part for 
the antisymmetric and symmetric components of the downwash velocity as 
determined from the Inner solution. Next, equations (13) are recombined to form 
corrected expressions for the source velocities on the airfoil and wake as 
follows : 


On the Ail. - oil: 


V^^corrected = V^, corrected + Vg, 

corrected 

(14a) 

V corrected = V , corrected * V , 

A o 

corrected 

(14b) 

In the Wake: 



corrected “ ^ 


(14c) 


A similar procedure is followed to correct the wake curvature term, 
leading to the expression 

dS = - tu^<D*+e) ^ (14d) 

where G„ is the ratio of the Inner solution to Its common part for the pressure 
lump across the wake streamline as determined from the trailing edge solution. 
Closed form expressions for the correction factors G , G , and G aa determined 
from the local trailing edge solution are given in reference 7. '^All three 
functions approach one and the matching conditions reduce to conventional 
boundary layer forms outside the trailing edge Interaction zone. Inside the 
Interaction zone equations (14) reduce to forms that will recover the local 
trailing edge solution when combined with solutions of the outer Invlscld 
equations. 

The correction functions completely eliminate the singularities arising 
In the conventional matching conditions and lead to solutions that are uniformly 
valid in the trailing edge region. We also note that the matching conditions 
lead to discontinuities in the pressure across the wake in the outer invlscld 
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a limited region near shock waves. Experience with the code indicates that 
the overall solution Is insensitive to the degree of local smoothing applied in 
this way. Although the theory is strictly rational only for airfoils with 
cusped trailing edges the program can be applied to airfoils with finite 
trailing angles. Good resi.lts have been obtained for airfoils with Included 
angles as larc^.e as 10°. 

The Iteration between the Inviscid and viscid solutions is controlled by 
monitoring the maximum residual, which is a parameter that measures the 
degree of convergence of the Inviscid solution. The boundary layer equations 
are solved and the matching conditions are updated whenever the maximum 
residual is reduced by a specified factor, typically equal to five. 

All calculations were carried out on an IBM 370/168 computer. Typical 
runs for a full viscous solution required about five minutes to reduce the 
maximum residual to 10*"^ on a 32 x 128 point grid. 

RESULTS 

In this section we report on results obtained with the new program. 
Comparisons of the theoretical results with wind tunnel data are presented 
for two airfoils, one a moderately rear loaded supercritical airfoil, the 
KORN I, which was tested at the National Aeronautical Establishment in 
Ottawa, and the other a newly designed NASA supercrltl:al airfoil with large 
rear loading that was tested in the NASA Langley 8~foo : transonic tunnel. 
Theoretical results are also presented which show the effect of neglecting 
the wake curvature term and of using a nonconservative formulation in the 
inviscid solution. 

The test of the KORN I airfoil at Ottawa was carried out at a Reynolds 
number of 21.7 x 10^ with the tunnel walls set at 20.5 percent porosity. 

Under these conditions blockage effects were small, but downwash effects 
due to the walls were large. Therefore, angle of attack corrections developed 
for the facility were applied to the data. (See reference 7 for further 
details.) The airfoil was aerodynamically smooth and was tested with natural 
transition. The calculations were carried out with transition fixed at 10% 
chord. This is a reasonable estimate considering the high Reynolds number 
of the test. The theoretical solutions were found to be relatively insen- 
sitive to the assumed position of the transition point. 

The theoretical surface pressures are compared with experrmental data in 
figure 4 for a free stream Mach number, M » 0.699 and an angle of attack 
a ■ 1.69° (corrected). The fully conservative solution and the experimental 
data are seen to be in good agreement. The agreement of the nonconservative 
solution with the data is much poorer. Since the drag coefficient of the 
nonconservative solution given in the figure does not Include the mass flow 
correction that Garabedian (ref. 13) has shown to be necessary, the favorable 
agreement of the uncorrected drag result with data is fortuitous. 
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correction of AM ■ -0.01 was employed In the theoretical computations. 

Although, accurate wall corrections are not available for this facility 
recent studies, carried out at NASA Langley Research Center, have indicated that 
corrections of this magnitude are appropriate. Transition was fixed in the 
tests by placing transition strips at 28 percent chord and was fixed at the same 
location in the computations. The tests were carried out at a Reynolds number 
of The results in figure 8 show that the new viscous theory 

with wake curvature not only predicts the shape of the drag polar but also the 
absolute levels of drag over a wide range of lift coefficients. The results 
also show that the neglect of wake curvature leads to a noticeable underpredic- 
tlon of the drag and that the nonconservative formulation leads to a signifi- 
cant overestimate of the drag. 

In figure 9 we compare the theoretical and experimental pressure 
distributions for a point about half way up the drag rise on the polar. The 
results show good agreement for both the pressure distribution and shock 
position. The slight underpredictions of the pressure over the rear upper 
surface were found to be due to numerical problems with code that have since 
been eliminated. 

The solution for the lift curve is given in figure 10. The solutions 
clearly show the large effect of wake curvature on the predicted lift 
coefficient. The type of differencing is seen to have only a slight influence 
on the predicted lift coefficient. Uncorrected data from the test are alsc 
plotted in the figure for reference. The results are suggestive of relatively 
large wall induced downwash on the airfoil. However since accurate angle 
of attack corrections were not available conq>arlsons with corrected data 
could not be provided. 


C(»)CLUSI0NS 

The main conclusions to be drawn from the results obtained to date with 
the new viscous theory are: 

o Wake curvature and trailing edge interaction effects are required 
in a completely consistent theory of viscous flow over airfoils. 

Together they account for about a half of the reduction of lift 
caused by the boundary layer. They are most important for rear 
loaded airfoils 

o The new theory ...eads to accurate predictions of lift and drag 
provided a fully conservative differencing scheme is used to 
solve the equations. Both the shape of the drag polar and the 
absolute levels of the drag are well predicted by the theory. 

The above are only tentative conclusions. The uncertainties in the drag 
measur'-ments of the Ottawa facility and the unknown level of the angle of 
attack and blockage corrections required in the Langley tunnel prevent us 
torn drawing more definite conclusions regarding the adequacy of the new theory. 

A more clear-cut validation of the theory would require improved wind tunnel 
tests with either small wall effects or accurate and documented corrections. 
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y = AU = -0^(6* + 0 )k 

Figure 1.- Flow structure and viscous matching conditions 



Figure 2.- Wake-curvature condition. 



a -.668 CO -.0062 CM --.0989 


Figure 3.- Outer inviscid solution for the pressure distribution whei 
CL denotes lift coefficient, CD denotes drag coefficient, and 
CM denotes pitching-moment coefficient. 



M - .699 ot » 1.69° RE - 21.5 X 10* 

CL - .668 CD - .0082 CM - -.0989 
CL - .654 CD - .0109 CM » -.1023 
CL - .669 CD - .0107 CM - -.1056 


Figure 4.- Pressure distribution on a KORN I airfoil where RE denotes 

the Reynolds number. 




Figure 8 



Figur- 7.- Drag polar for the KORN I airfoil. 
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• - Drag polar for the NASA supercr Ideal airfoil 
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